Interplay between quantum dissipation and an in-plane magnetic field in the spin ratchet effect 
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We investigate the existence of the pure spin ratchet effect in a dissipative quasi-one-dimensional system with 
Rashba spin-orbit interaction. The system is additionally placed into a transverse uniform stationary in-plane 
magnetic field. It is shown that the effect exists at low temperatures and pure spin currents can be generated by 
applying an unbiased ac driving to the system. An analytical expression for the ratchet spin current is derived. 
From this expression it follows that the spin ratchet effect appears as a result of the simultaneous presence of 
the spin-orbit interaction, coupling between the orbital degrees of freedom and spatial asymmetry. In this paper 
we consider the case of a broken spatial symmetry by virtue of asymmetric periodic potentials. It turns out that 
an external magnetic field does not have any impact on the existence of the spin ratchet effect, but it influences 
its efficiency enhancing or reducing the magnitude of the spin current. 
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I. INTRODUCTION 

It is well known that a directed stationary flow of particles 
in a system can be created by unbiased external forces. In 
general this possibility arises when the system is not invariant 
under reflections in real space. This fact is mainly indepen- 
dent of the mechanics which underpins the particle motion, 
classical or quantum. However, the microscopic origin of this 
effect, conventionally called the ratchet effect, is different in 
the classical and quantum case. One principle source of that 
difference is quantum mechanical tunnelling which does not 
have analogs in the classical mechanics. Correspondingly, one 
usually distinguishes between classical and quantum ratchet 
effects. In this paper we concentrate on the latter one in a dis- 
sipative system. Such dissipative ratchet systems act as Brow- 
nian motors 12 turning Brownian into directed motion. The 
existence of the ratchet effect in a quantum dissipative one- 
dimensional (ID) system which lacks the spatial symmetry 
has been first theoretically predicted in Ref.[3l Later, within a 
tight-binding model where the lowest bands are narrow, it has 
been disclosed that a ratchet state of the particle transport can 
only be achieved when at least the two lowest Bloch bands 
contribute to transport 4 . To obtain the ratchet effect in sys- 
tems with weak periodic potentials at least two harmonics of 
the potential should enter the dynamical equations 5 . Rectifica- 
tion can also take place in a single-band tight-binding model 
where the spatial asymmetry is concealed from the electron 
dynamics. One way to achieve this is to use unbiased external 
forces with harmonic mixing 6 . 

Coherent charge ratchets based on molecular wires with an 
asymmetric level structure of the orbital energies were pro- 
posed in Ref. 7. In this case weak dissipation originates from a 
weak coupling between the wire edges and leads which repre- 
sent fermionic reservoirs. In contrast to the systems described 
above in this system there is no dissipation in the wire. The 
ratchet effect is a result of the dissipative coupling of the wire 
to fermionic baths. 

In a different branch of condensed matter a new research 
field has emerged during the last decade, namely spintronics, 
where one tries to make use of the spin degree of freedom of 



a particle instead of only the charge one. One essential differ- 
ence between spin and charge is that a particle can have more 
than one spin state while it has only one charge state. In the 
context of transport it is important that the spin state of a parti- 
cle can strongly depend on the transport conditions, in particu- 
lar on the transport direction, as it happens for example in sys- 
tems with spin-orbit interaction. This fact has founded a new 
arena for different spin devices used to store, transform and 
transfer miscellaneous information. The possibility to transfer 
the spin separately from charge plays an important role. This 
can be implemented by so-called pure spin currents, that is 
spin currents which are not accompanied by charge currents. 
Thus the generation of such currents has been extensively 
discussed. Among different mechanisms of spin-orbit inter- 
action Rashba spin-orbit interaction (RSOI) 8 plays a distin- 
guished role because it provides an opportunity to control the 
spin-orbit coupling strength by an external electric field. The 
change in the band structure spawned by the spin-orbit inter- 
action leads to one of the most remarkable effects in spintron- 
ics, the intrinsic spin-Hall effect, first predicted by Murakami 
et al~ for hole-doped semiconductors with the spin-orbit in- 
teraction of the effective Luttinger model for holes and later 
by Sinova et al— in a high-mobility two-dimensional elec- 
tron gas (2DEG) with RSOI. The spin current which results 
from the intrinsic spin-Hall effect is pure and its experimental 
detection was discussed, e.g., by Wunderlich et al.- Another 
kind of spin-Hall effect, the extrinsic spin-Hall effect, is a re- 
sult of the spin-orbit interaction as well. The spin currents 
related to the extrinsic spin-Hall effect are also pure. Such 
pure spin currents were experimentally detected through op- 
tical measuring of electron spin accumulation at the edges of 
the samples^ and through the reciprocal spin-Hall effec t 13 ' 14 
in Ref. 15. Another approach to create pure spin currents is to 
use polarized light. For example in noncentrosymmetric semi- 
conductors one-photon absorption of linearly polarized light 
induces pure spin currents 16 . The pure spin current response 
to linearly and circularly polarized light irradiation, exciting 
electrons from valence bands into the conduction bands, was 
studied by Li et al— and by Zhou et al. 18 for 2DEGs with 
RSOI. An alternative technique of getting pure spin currents 
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is quantum spin pumping. The idea of quantum spin pumping 
comes from the general idea of electron pumping 19 . Electron 
pumping assumes that in a given system any voltage bias is 
absent and the particle flow is a result of a cyclic variation of 
at least two parameters of this system. When the electron spin 
is involved due to some mechanisms, various quantum spin 
pumps emerge. For example spin pumps based on electronic 
interactions 20 , magnetic barriers 21 , carbon nanotubes 22 have 
been discussed. A spin pump based on a quantum dot was 
experimentally implemented by Watson et alP The pure spin 
current generation using the spin ratchet effect in coherent and 
dissipative systems with RSOI was investigated in Refs. [24] 
and|25| respectively. The spin ratchet effect in the presence of 
a non-uniform static magnetic field without spin-orbit inter- 
action, the Zeeman ratchet effect, was studied in Ref. [26| for 
coherent quantum wires formed in a 2DEG. However, the spin 
ratchet effect in a dissipative system in an external magnetic 
field has not been considered up to now. 

In this paper we extend the results of Ref. [25| to include a 
transverse in-plane uniform stationary magnetic field. Specifi- 
cally, we consider non-interacting electrons in a quantum wire 
formed by a harmonic transverse confinement in a 2DEG with 
RSOI. The electrons are also subject to a ID periodic poten- 
tial along the wire direction and the in-plane magnetic field 
perpendicular to the wire. An orbital coupling between this 
originally isolated system and an external environment causes 
dissipative processes affecting indirectly the spin dynamics 
through RSOI. 

An external ac driving originates in our work from an ap- 
plied ac electric field. We show that for such a driving the net 
stationary charge current is strongly suppressed if the trans- 
port is governed only by electrons of the Bloch sub-bands re- 
lated to the same Bloch band which would result from the 
corresponding truly ID problem without RSOI. However, at 
the same time and under the same conditions a net stationary 
spin current turns out to be activated in a spatially asymmet- 
ric situation and for finite values of the spin-orbit coupling 
strength and the coupling strength between the orbital degrees 
of freedom. The magnetic field does not destroy this picture, 
but it can partly reduce or on the contrary enhance the ratchet 
effect. 

The paper is organized as follows. In Section|II]we describe 
a model within which a ratchet like behavior of the spin trans- 
port can be achieved and present a master equation in terms of 
populations and transition rates between the basis states used 
to calculate the charge and spin currents. These basis states 
are then thoroughly discussed in Section [TTTJ A tight-binding 
model is formulated in Section [IV] In Section [V] we present 
the transition rates and their properties. Finally, in Section 
IVTl we derive analytical expressions for the charge and spin 
currents and explore the spin ratchet effect in the system. 



H. FORMULATION OF THE PROBLEM 
The full Hamiltonian of our problem is 

fffullW = H + H ext (t) + -ffbath, (1) 



where H is the Hamiltonian of the isolated periodic system, 
H ext (t) describes an external driving and Hbnth represents the 
term responsible for dissipative processes. 

The isolated quasi- ID periodic system is formed in a 2DEG 
(x — z plane) with RSOI using a periodic potential along the x- 
axis and a harmonic confinement along the z-axis. The whole 
system is in a uniform stationary magnetic field along the z- 
axis: 
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where Hq is the z-component of the magnetic field Ho = 
(0, 0, Hq), and we have used the gauge in which the compo- 
nents of the vector potential are A x — —Hoy, A y = A z = 
(Landau gauge). Additionally, we have taken into account the 
fact that in a 2DEG y = 0. In Eq. (O the operator k is related 
to the momentum operator p as p = fik, qjq is the harmonic 
confinement strength, k so the spin-orbit interaction strength, 
7 the strength of the coupling between the orbital degrees of 
freedom x and z, g the electron spin g-factor, /.ib the Bohr 
magneton, and U(x) denotes the periodic potential with pe- 
riod L, 



U(x + L) = U(x). 



(3) 



In the following we assume that the periodic structure is 
subject to an external homogeneous time-dependent electric 
field. Only the x-component of the electric field vector is non- 
zero, that is the electric field is parallel or anti-parallel to the 
x-axis. Experimentally this can be implemented using for ex- 
ample linearly polarized light. The external force thus couples 
only to the x-component of the electron coordinate operator: 



H, 



-F(t)x, 



(4) 



where the force F(t) is unbiased. In this work we use the 
time-dependence 



F(t) = Fcos(n(t-t )). 



(5) 



The term "unbiased external force" should not be confused 
with voltage bias. An external force is called unbiased if it 
is periodic in time and its mean value, that is its average over 
one period, is equal to zero. It is obviously our case as one 
can see from Eq. (0. 

The system is also coupled to an external bath. In the 
present work we assume the transverse confinement to be 
strong enough so that the probabilities of direct bath-excited 
transitions between the transverse modes are negligibly small. 
In other words, the wire is truly ID from the point of view 
of the bath which directly changes only the dynamics along 
the wire. Thus in our model the external environment cou- 
ples to the electronic degrees of freedom only through x. The 
bath itself as well as its interaction with the quantum wire are 
described within the Caldeira-Leggett mode l 27 ' 28 , 
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The bath is fully characterized by its spectral density defined 

as 

No 2 

j{u) = !r-5_f( u - Ua ). (?) 

a— 1 

It is important to emphasize that, due to the spin-orbit interac- 
tion and orbit-orbit coupling, the direct dissipative interaction 
between the longitudinal dynamics in the wire and the exter- 
nal environment has an indirect impact on the transition rates 
between different transverse modes. The transverse dynamics 
in the wire indirectly feels the presence of the external bath 
through the spin-orbit interaction and orbit-orbit coupling. 

The dynamical quantities of interest are the charge and spin 
currents. Specifically, the longitudinal charge current Jc{t) 
is given (see for example Ref. 4) as a statistical average of 
the longitudinal charge current operator Jc(t), i.e. the prod- 
uct of the velocity operator v(t) and the elementary electronic 
charge — e, 

J c (t) = -efi(t), (8) 

Mt) = -ejTv[xp(t% (9) 

where p(t) — Trb a thW / (0 is the reduced statistical operator 
of the system, that is the full one W(t) with the bath degrees 
of freedom traced out. 

For the longitudinal spin current operator we use the defi- 
nition suggested by Shi et al.—, 

Jh{t) = j t {a i x), (10) 

which was further developed and applied to a two- 
dimensional hole gas by Zhang et al^- The advantage of this 
definition over the conventional one ( Jg = o^t)) is that using 
the corresponding spin current, 

4(t) = fax p(t)), (ID 

the continuity equation for the spin density can often be writ- 
ten without a source term, which means that the spin current 
defined in this way is conserved. This conserved spin current 
can be uniquely related to the spin accumulation at a sample 
boundary. The out-of -plane polarized spin accumulation can 
experimentally be measured with Kerr rotation microscopy 31 
or the Faraday rotation technique 32 . The in-plane spin polar- 
ization is not directly measured by Kerr rotation microscopy, 
but it can still be scanned by a magneto-optic Kerr microscope 
using, e.g., the cleaved edge technology as discussed by Ko- 
tissek et aZ. 33 Even when the continuity equation contains a 
source term, there is still one advantage of the spin current 
operator definition (TlOb . This definition leads to a very rea- 
sonable physical result: the corresponding spin current in (fTTT i 
vanishes in insulators. In Section [VI] we will return to this 
point and analytically prove that when the periodic potential 
gets stronger and as a result the energy bands get narrower, 



that is when the system turns into an insulator, the spin cur- 
rent given by Eq. (TTTb goes to zero. Below we will calculate 
only the spin current polarized along the z-axis and denote 
this current as Jg, i.e., Js{t) = Ji(t). The components of 
the spin current polarized along the x and y axes are zero as 
shown in Appendix [B] The discussion of the difference be- 
tween the conventional definition of the spin current and the 
spin current definition used in our work can also be found in 
Appendix 151 

It is convenient to calculate the traces in (0 and (TTTb us- 
ing the basis which diagonalizes both x and a z , because this 
requires to determine only the diagonal elements of the re- 
duced density matrix. In a quasi- ID periodic system with 
RSOI the energy spectrum can be related to the one of the 
corresponding truly ID problem without RSOI 34 . This links 
the Bloch bands of that truly ID problem to the Bloch sub- 
bands of the quasi- ID problem. The general structure of the 
results obtained in Ref. [34] is retained in the presence of the 
orbit-orbit coupling and a uniform stationary magnetic field 
along the z-axis. A slight change of the theory is given in 
Appendix ICl We shall consider electron transport under such 
conditions when only a finite number of the Bloch sub-bands 
is involved. The basis which diagonalizes x and a z becomes 
in this case discrete. The total number of the Bloch sub-bands 
is equal to the product of the number, Nb, of the Bloch bands 
from the corresponding truly ID problem without magnetic 
field and without spin-orbit coupling, the number, N t , of the 
transverse modes and the number of the spin states. In this 
work we shall use the model with N-q — 1, Nt — 2. Since 
there are only two spin states, the total number of the Bloch 
sub-bands in our problem is equal to four. The representation 
in terms of the eigen-states of the coordinate operator for a 
model with discrete x-values is called discrete variable rep- 
resentation (DVR)2£. Let us call c-DVR the representation 
in which both the coordinate and spin operators are diagonal. 
Denoting the a-D VR basis states as { | a) } and eigen-values of 
x and a z in a state |a) through x a and a a , respectively, the 
charge and spin currents (O and ( TTTb are rewritten as 

Jc(t) = -e^2x a -^P a (t), 

" d (12) 

■^s(t) = y^^ a x a —p a {t), 

a 

where P a (t) = (a\p(t)\a) is the population of the ct-DVR 
state \a) at time t. 

We are interested in the long-time limit of the currents 
Jc{t) an d Js(t) averaged over the driving period T = 2tt/Q 
with the time average of a time dependent function f(t) de- 
fined as /(<) = (1/T) J t t+T dt'f{t'). From (EJ) it follows 

Jc(t) = -ey^Xq-^-Pqft), 

d - (B) 
Js(t) = X! o- a x a —P a (t). 

a 

The advantage of working in the cr-DVR basis is that real- 
time path integral techniques can be used to exactly trace out 
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the bath degrees of freedom 36 ' 37 . Moreover, at driving fre- 
quencies larger than the ones characterizing the internal dy- 
namics of the quasi- ID system coupled to the bath, the aver- 
aged populations P a (t) can be found from the master equa- 
tion, 



dt 



p a (t)= ivp/jm- f 0ocP a {t), d4) 



P 

(P^c 



valid at long times. In Eq. ( TBI T a p is an averaged transi- 
tion rate from the state \j3) to the state \a). In order to obtain 
concrete expressions for the averaged currents the rr-DVR ba- 
sis must be found explicitly. This is the subject of the next 
section. 



III. DIAGONALIZATION OF a z AND x: THE <r-DVR 
BASIS 

The eigen-states of the a z operator were found in Ref. [34] 
(see Eq. (12) therein) for a model without coupling be- 
tween the orbital degrees of freedom and magnetic field. The 
changes necessary to include those two effects are discussed 
in Appendix|C] The eigen-value equation for the a z operator 
is 



a z \l,k B ,j,cr) 7 ,j = a\l,k B ,j,cr)-y t j. 



(15) 



In Eq. (15[ I, k B , j, u stand for the Bloch band index, quasi- 
momentum, transverse mode index and z-projection of the 
spin, respectively. Since in the presence of the orbit-orbit cou- 
pling the periodic potential U 7 j(x) (see AppendixO depends 
on 7 and j, we have labeled the ket-symbol with the subscript 
7, j. In the ensuing analysis we follow the same rule and la- 
bel all the bra- and ket-symbols with the subscript 7, j, that is 
7J (---|and|---) 7J . 

It is convenient to start the diagonalization of the coordinate 
operator writing its matrix in the {\l, k B ,j, cr) 7 .j} representa- 
tion: 



7,j'C i^Bii j a 2-1') ^Bj j, <7 )-/,j — 

The diagonal blocks, 

7,iv ! '> fe B>.7= cr = l|£|^ fc B, j,<7 = l) 7 ,j = 

= (J, , k B + fc so |a;|/, k B + k BO )j t j, Vj, 

7 ,j i ^Bi h a = ~~ 1 l^l'i ^B: j, a = ~~ 1)7 j — 
— , fcg — k so \x\l,k B — fc so ) 7 j, Vjj 



(16) 



(17) 



are unitary equivalent for a given value of the index j and thus 
the eigen-values of x do not depend on a. 

As it is shown in Appendix lAl the eigen-values of the ma- 
trix k' B \x\l, k B )j,j are 



(18) 



where m = 0, ±1, ±2 . . ., £ = 1,2,..., N B and the eigen- 
values cL ; fj are distributed within one elementary cell. If, 



for example, the system is divided into the elementary cells 
in such a way that the origin of coordinates is at the center 
of an elementary cell, then —L/2 < d 7; £j L/2. In Eq. 
(fT8l we have taken into account that the periodic potential 
U 7 j(x), introduced in Appendix ICl depends on 7 and j, and 
thus the eigen-values distributed within one elementary cell 
also acquire a dependence on 7 and j. 

From (TToT l and ( [T8"b it follows that one can label the eigen- 
states of x with the quantum numbers £, m, j, a, that is 
as |£, m,j, er) 7 j, and in the {\l, k B , j, c) 7 .j} representation 
these eigen-states have the form: 



7,J' (*) fc B, /, cr'IC, m, j, cr) 7 j = 



(19) 



The corresponding eigen-values are x 7; <7 m ,j iC7 = x^-^^mj. 
From the eigen-value equation 



(20) 



written in the {\l,k B ,j, c) 7 j } representation through the use 
of©, 



^ 7J -(i, fc B + crfc so I i 1 1 ! 



crfc s . 



(21) 



(22) 



X 7j( ? : fc B 5 J> Cr IC,^,i,C r ) 7 j = 

= x 7;C,m ,j 7J'('' fc B,j, cr|C,m, j, er) 7J -, 
it follows that 

7,j(^ fc B, j, CT = l|C,m,j,(r= 1) 7J - = 

= k B + fc so |Cl m )l,jl 

~,,j{l,k B J,CT = -l|C,m,i,CT = -l) 7 ,j = 
= jj(l 7 k B fc so |£, TTl)jj. 

Since |£, to, j, er) 7 j is also the eigen-state of <r z corre- 
sponding to the eigen-value a^ jm> j t(T = a, we infer that the 
cr-DVR basis states |a) from the previous section are just the 

|C,m, j,(r) 7 ,j states, that is {|a}} = {|C, m, j, <r) 7J -}. 



IV. (j-DVR REPRESENTATION AND ITS TIGHT-BINDING 
MODEL 

Let us represent the Hamiltonian H in the er-DVR basis 
obtained in the previous section in order to derive an effective 
tight-binding model. 

Using the { |£, m, j, cr) 7 ,j } basis the Hamiltonian H can be 
written as 

H= ^2 ~f,j'(C m ', f ,v'\H\C,m,j,(j) ltj y. 



C',m',/V 



(23) 



x |C ,m ,j ,0 } 7i f Jtj ((,m,j,a\, 
with the matrix 
7d'(C> , m'j' ) a'\H\C ) m,j,a)y tj = e 7; (^(fc B )x 

x j,f(C,m'\l,kB +o-'fc so ) 7 j/X (24) 

(i',^')x 



5 



The tight-binding approximation of ( l23T > is obtained if one we finally have 



assumes that the matrix elements (1241 with \mf — m\ > 1 are 
negligibly small. 

We consider temperatures low enough and assume that 
electrons populate only the lowest Bloch sub-bands with / = 1 
(i.e., Nb — 1)- Under this condition the periodic potential can 
be of arbitrary shape and the only limitation on it is the valid- 
ity of the tight-binding approximation. 

Below we thoroughly study the case where the four low- 
est Bloch sub-bands are the ones with / = 1, ry = 1,2,3,4 
and the only ones which are populated with electrons. For 
simplicity we consider weak orbit-orbit coupling and calcu- 
late the corresponding eigen-energies e^^fes) and eigen- 
spinors T ,i.k B ,ri(j: f) retaining only the first two transverse 
modes, that is j — 0, 1. In this case H has the form 

H = ^2 ^£ 7 y )CT |m,j,o-) 7> j l!j (m,j,a\+ 



3, 

■ \ ^ a intra 



t(r (m)\m, j',(r') 1 ,j' 7 ,j(m, j,o\ + 



(25) 



where 



|w,i,(r) 7 j = |C= l,m, j,a) r 



(26) 



and we have defined the on-site energies e 7; j,cr and hop- 
ping matrix elements A^^ a ,. ja (m), A™^,.^ a (to) and 
A™*f£., ff (tn) as follows ' 



7;i'i°'';ji | 



(m, j, cr|i?|m, j, <t). 



A^V ;i , CT (m) = 7J v(m,j',CT , |i?|m,i, ( j) 7J , 

A inter, b / \ 
., , to = 



7U ,<r \3-.' 
, inter, f 



A 1 ', (ml 



*1\3 ,<* \3,' 

Note that 



■ (m,j',a'\H\m + l,j,a) l!j , 
>(m+ l,j',a'\H\m,j, cr) 7J -. 



(27) 



T,3 



(to), 



r a inter, b / \-\* 

Introducing the notations 



A inter,f , s 



(28) 



(29) 



{O = {(?»}, 

£ = 1^(0,1), £ = 2^(0,-1), 
^ = 3^(1,1)^ = 4^(1,-1), 



(30) 



•£=i 

4 

X Ai;f; a e (™)|™,^) 7 ,«' 7 ,«(m,ei+ 

^«'=i (31) 
£ (A^ b (m)|m,f') 7 ,e 7 ,«(m + l,e| + 

£,£'=1 V 



+ A^f/(m)|m + l,0 7 , e 7 , e (m,£| 



Equation d3TT > represents a tight-binding model which can now 
be used to perform actual calculations of quantum transport in 
a dissipative system. 

To conclude this section, we would like to note that because 
of the simultaneous presence of the harmonic confinement and 
RSOI the system splits into two subsystems. The first sub- 
system is characterized by £ = 1,4 and the second one by 
£ = 2,3. These subsystems are totally decoupled: there is no 
electron exchange between them. Such a state of affairs per- 
sists if one considers more than two transverse modes. In this 
work, for simplicity, we only consider one subsystem, namely 
the one with £ = 1,4. Such uncoupled subsystems also appear 
within the hard wall model of the transverse confinement^. 



V. TRANSITION RATES 

The tight-binding model introduced in Section [IV] relies 
upon the fact that the hopping matrix elements (l27i i are small. 
In this case the second-order approximation for the averaged 
transition rates in Eq. ([Pil l can be used givingi^i 



I a m ,m 1 2 



A 7;£',€ _ 
dre -[(x 7 . m , ? -x„ 



h 2 



% , i(r ) 2 mQ(T)+i[(e y . i -e Ti(r )/H]T x (32) 



X Jn 



,e) ■ (fa 

sm( — 



where x T . m ^ = x 7; £=i, m ,£ = mL + d 7; £ with d 7; £ = d 7; i j, 
A™^',™ = 7l £/ (to', £'\H \m, £) 7 ,£ the hopping matrix element 
between the states |to',£') 7 .,c/ and \m, £) 7 ,f, Jo(%) the zero- 
order Bessel function and Q(t) the twice integrated bath cor- 
relation function^: 



1 r°° 

Q(t) = - / dui 
k Jo 



J(w) 



x [1 — cos(wr)] + isin.(wr) 



(33) 



where J(ui) is given by Eq. (0 and j3 is the inverse tempera- 
ture. 

The transition rates are functions of the orbit-orbit cou- 
pling strength 7 because the Bloch amplitudes as well as 
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the difference Ad 7 = g? 7; i.o — <^ 7 ;i,i depend on 7. Within 
the context of the tight-binding model the eigen-values d 7; i.o 
and d Tt i t % tend to zero and fulfil Ad 7 /Z r <C 1, where l r = 
min[X, ^/S/muo, Ml/F, . . .]. Consequently, the transition 
rates depend on 7 predominantly through the Bloch ampli- 
tudes, and in this work we pay no regard to terms of order 
0(Ad J /l I ). This is also consistent with our model taking into 
account only the first two transverse modes. Keeping terms 
of order 0(Ad 7 /^ r ) would mean that the strength 7 of the 
orbit-orbit coupling is large enough so that one would need to 
consider more than just the first two transverse modes because 
in this case the non-diagonal elements would be comparable 
with the diagonal ones. 
Using the notations, 



-\m,m pmtra 



=\m,m-\-l 



pintcr,b 

:i 7 ;£',« 



eve, 



(34) 



p?n+l,m T-dnter,f 



from d32t one obtains 



A™^ r t b (m) and A™^ f (to) only through a phase factor as 
it is shown in AppendixlDl From ( 1291 and d36l i it follows that 



pintcr,b pintcr,f 



pintcr,bpintcr,b pinter,f pintcr,f 



(38) 



(39) 



To calculate the charge and spin currents we additionally 
need the transition rates 

P pintcr, f . pintra , fainter, b /Al\\ 

1 7 ;W = 1 T ,w + 1 7 ;S,«' + 1 T ,t,e ■ (40) 

As pointed out at the end of Section IIVI the system is split 
into two subsystems isolated from each other. Since electron 
exchange between the subsystems is absent one can write 



r 7; i,2 — r 7; i,3 — r 7; 2,i — r 7; 2,4 — 
= r 7: 3 1 = r~ ; 3 4 = r 7: 4 2 = r~ ; 4 3 = 



(41) 



and 



pintra 



(35) 



The last equalities are very useful because they allow us to 
significantly simplify the expressions for the charge and spin 
currents which are derived in the next section. 



where 

x J 



inter, b 

r,6',S 

inter, f 

r,6',(. z 



:|A^ b (TO)|V 7;e ,e, 

|A^ f (m)| 2 J 7;e , £ , 



(36) 



h 2 



2FL ( fir 



(37) 



Note that ri^f^ b and T 1 ^ 1 '} do not depend on to due 
to the Bloch theorem which leads to an rn-dependence of 



VI. CHARGE AND SPIN CURRENTS 

The expressions for the stationary averaged charge and spin 
currents, 



J, 



c = lim J c (t), J s °° = lim J s (t), 



(42) 



can be found from the averaged master equation (fT~4T > which 
we rewrite here using the er-DVR indices and tight-binding 
approximation introduced in Section [IV] and utilizing the no- 
tations of Section[V]for the transition rates: 



,4 4 

r>m /.\ \ * [pinter,f r>m-l i ,\ . pintra r>m / + \ , pinter.b D m+1 / , \] \ "> [pintcr, b . pintra , pintcr, f] rjm 



?'=i «'=i (43) 

winter, f 73m — 1 /«\ , fainter, b nm+1 /ja] [fainter, b . pintcr, f 



pmter,f j^m — 1 / ,\ . pinter,b r,m+l / ,\1 fpinter,b . pintcr, f] j^m />\ 



From (fT2l and (l43l One finds Too r/j j \ /pinter.f , winter, b\ , 

, 7- /fainter, f fainter. b\ 1 oo 



4 where we have used Eq. ( fT8l . To derive Eq. d45t we have 

= -eh [r^e - f 3$ b ]p~ fi / , (44) additionally made use of Eq. g5J. In Eq. g5) = «7 C=W 
f g/=i and o"i = (T3 = 1, (72 = ^4 = — 1 as it follows fromEq. ( f30b . 
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The quantities are defined as 

p 7;f (t) = V^(4 p£ f = limp 7i £(i), (46) 

m 

and they satisfy the constraint 

p T ,l(t) + p Ti2 (t) + p Tt3 (t) + p rA (t) = 1, Vt. (47) 

As already mentioned at the end of Section |IV] we only con- 
sider the subsystem with £ = 1,4. The properties of the sta- 
tionary averaged transport do not depend on initial conditions. 
We choose the following ones: 



p 7;1 (i = 0) = l, p TA (t = 0)=0. 



(48) 



Because of the constraint ( |47T > p T ,2(t = 0) = p T .3(t = 0) = 
and since there is no electron exchange between the subsys- 
tems, the states of the subsystem with £ = 2, 3 remain empty 
at any time, p-y-^it) — Pj-3(t) = 0, Vi. This leads to 
P^2 — Pt°3 = 0- Then from the master equation d43l > with 
the initial conditions d48b and using (|40| >. (|4T) one obtains 



Pj;l 



TAA 



T 7; l 



tAA 



- 7;4,1 



P-y-A 



r 



t;4,i 



r 7; i,4 + r 7; 4 ; i 



(49) 



Using Eqs. ([35}, (OS"© and (|49) it follows from tO 



0, 



(50) 



that is the absence of the stationary averaged charge transport. 
However, using Eqs. ((35), d38) . d40l > and (|49) we get from Eq. 



TOO 

J s — 



2L 



r 7; i,4 + r 7; 4,i 



/pintcr,f pintcr,b T^intcr,bTnintcr,f \ /ci\ 
l i 7;l,4 i 7;4,l — 1 7;!, 4 1 7; 4,1 J- ^U 



The last expression can be rewritten in terms of the hopping 
matrix elements A™!^ (m). Making use of Eqs. ((29), ((35), 
and (|40) we derive the stationary averaged spin current: 



joo _ ^lAA'^iAA 



X 



J~f;lA + JfAA 

x(|A^ f (m)| 2 -|A^ r /(m)| 2 ) 



(52) 



Using Eqs. d52) and (1D21 > the non-equilibrium stationary av- 
eraged spin current can be written as 



'n-e,S 



T T 7 (0) 7 (0) 

Jf-,1, 4^7:4,1 J 7;l,4' y 7;4,l 



77? 



J 7 .! 4 + J 7 . 4 ! t(0) , t(0) 
7,1,4 7,^,1 J 7; i^4 i" J 7; 4,l 

sin[(fc B - fcg)L] Im[F 7;fcBifc; 



(53) 



where j' l f is given by Eq. ((37) with F = and the func- 



tion F 7; fc B) fc^ is defined by Eq. (1D31 >. Note the structure of Eq. 
((53). It is the product of two factors of different physical ori- 
gin. The factor in the second line describes the isolated system 
and the factor in the first line is purely due to the dissipative 



coupling to an external environment. To get Eq. (153) we have 
eliminated from the equilibrium spin current arising due 
to the non-compensation 40 of the spin currents from different 
bands of the Rashba-Bloch spectrum of the isolated system. It 
turns out that this effect is strong enough to indenture in a dis- 
sipative system. Below we only consider the non-equilibrium 
spin current, J^L C s , and not the full one, J^°. 

Let us at this point also mention the dependence of the 
spin current J^L C s on the magnetic field Hq- Since the mag- 
netic field is applied along the z-axis, it couples to the system 
through the a z operator and thus the hopping matrix elements 
y. (m) do not depend on Hq. It then follows that the 



inter, f 
7;1(4),4( 



spin current depends on the magnetic field only through its 
dissipative prefactor. The dependence on Hq comes into play 
through the on-site energies e 7; i(4). The difference e 7; 4 — e 7; i 

which enters the integrals J 7; i(4),4(i) and J^lu^ 4m) 
written as: 



can be 



-7;4 



-7;i 



k E 

hujQ 



(54) 



2g^ B H 



where N is the number of the elementary cells and e^.^fcs) 
are the eigen-values of the truly ID Hamiltonian 



ID _ 

0;7J — 



h 2 k 2 x 

2m 



U(x) 



'7 



h f . 1 

mujoL 2 \ 2 



(55) 



Therefore, in the presence of a transverse in-plane uniform 
stationary magnetic field the existence of the spin current is 
possible under the same conditions which were discussed in 
Ref. 25[ For completeness we list these conditions below. 

From d53) one finds, as mentioned in Section [TTJ that when 
the electronic states become localized, the stationary averaged 
spin current vanishes. Indeed, in this insulating limit the func- 
tion F^kntf does not depend on the quasi-momenta fee and 
k' B and Eq. ((53) gives zero. 

When the spin-orbit interaction is absent, that is fc so = 0, 
we get from I 



n-e,S |fc 



,= 0. 



(56) 



Further, if the orbital degrees of freedom x and z are not 
coupled, that is 7 = 0, it follows from Eqs. ( 153) and (ID6t that 



J, 



n-c,S 



,= 0. 



(57) 



Finally, if the periodic potential is symmetric, the Bloch 
amplitudes are real and we find from Eqs. (153) and ( ID4I ) 



J, 



n— e,S 



= 0, for symmetric periodic potentials . (58) 



Summarizing the results of this section we conclude that in 
order to generate a finite stationary averaged spin current three 
conditions must simultaneously be fulfilled: 1) presence of the 
spin-orbit interaction in the isolated system; 2) finite coupling 
between the orbital degrees of freedom x and z; 3) absence of 
the real space inversion center in the isolated system. 
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Among these three conditions the second condition is per- 
haps less transparent and a simplified physical interpretation 
is necessary. We propose the following physical explanation. 
The orbit-orbit coupling leads to the situation in which the 
strength of the periodic potential varies across the quasi- ID 
wire. The periodic potential is equal to U (x) in the center of 
the wire and gets stronger closer to its edges. Thus the elec- 
tron group velocity is larger in the central region of the wire 
and decreases closer to the edges. At the same time the elec- 
tron distribution across the channel depends on the transverse 
mode j. It is given by the Hermite polynomials. For j = the 
electrons populate the center of the wire while for j = 1 they 
are distributed in regions closer to the edges. Hence, the elec- 
trons with j = are faster than those with j = 1. Because of 
the mixing between the confinement and RSOI different trans- 
verse modes carry different spin states. Therefore, we con- 
clude that different spin states have different group velocities 
along the wire. This difference results in a finite longitudinal 
spin current. 

Finally, one observes that a transverse in-plane uniform sta- 
tionary magnetic field alone is not enough to produce the spin 
current in a driven dissipative system. The magnetic field can 
only affect the magnitude of the spin current when the prop- 
erties of the isolated system meet the three conditions derived 
above. 



VII. RESULTS 

In this section we show some results obtained numerically 
for the theoretical model developed in the previous sections. 
As an example we consider an InGaAs/InP quantum wire 
structure. The values of the corresponding parameters used 
to get the results are similar to the ones from the work of 
Schapers et al£^ In particular, Hujq = 0.225 meV, a = 
h 2 k so /m = 9.94xlCr 12 eV-m (which gives fc so = 4.82 xlO 6 
m _1 ), rn = 0.037mo (mo is the free-electron mass). The 
value, g — 7.5, of the electron spin g-factor (in our nota- 
tions g = —g*/2, where g* is the effective gyroscopic factor 
measured experimentally) is taken from Ref. 42. From these 
parameters and for example for the period of the super-lattice 
L = 2.5^/h/mujo « 0.24 /im, which is easily achievable 
technologically at present 43 , it follows that k so L « 0.3687T. 



The asymmetric periodic potential is 



U(x) = hiuo< 2.6 



1 



cos 



(2ttx 



1.9 



1.9cos( — ) 



(59) 



The bath is assumed to be Ohmic with exponential cutoff: 

J(uj) — T^expf — — ) , (60) 

where r\ is the viscosity coefficient and lo c the cutoff fre- 
quency. We use uj c = 10 ojq. 
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Amplitude of the driving force, F [ h d>JL ] 



FIG. 1: (Color online) Non-equilibrium spin current, J^L e s , as a 
function of the amplitude, F, of the driving force for different val- 
ues of the z-projection of the magnetic field Hq. Further parameters: 
temperature fceoitz.r = 0.5, spin-orbit coupling strength k so with 
k so L — 7r/2, orbit-orbit coupling strength 7 = 0.08, driving fre- 
quency Q = 0.2, viscosity coefficient 77 = 0.08. 
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Spin-orbit coupling strength, k^LIn 

FIG. 2: (Color online) Non-equilibrium spin current, J^L e ,S' as a 
function of the spin-orbit coupling strength, fc so , for different values 
of the ^-projection of the magnetic field Ho- The driving amplitude 
is F — 1.0 hujo/L. The other parameters are as in Fig.Q] 



To present the results we use in all the figures the units of 
hujo and ojq for energies and frequencies, respectively. The 
viscosity coefficient is taken in units of txiloq. 

Let us discuss possible values of the driving parameters. In 
a dissipationless system (or in a system with weak dissipation) 
of size L one should restrict possible values of the driving 
amplitude and frequency, < FL < Hujq and < £1 < ojo, in 
order to stay within the validity of the model with the first two 
transverse modes opened. In a strongly dissipative system, as 
in our case, it is not necessary to fulfil the last inequalities 
because an electron loses a huge amount of its energy due 
to intensive dissipative processes. In general, our model of 
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a driven strongly dissipative system taking into account the 
first four Bloch sub-bands remains valid if at long times the 
electron energy averaged over one period of the driving force, 
e av (F, ft, rj) (which is a function of the driving and dissipation 
parameters), is smaller than ftioo, e a v(F,£l,r]) < Huq. This 
can take place even if FL > ftwo and > ujq because even 
at such driving the strong dissipation (large values of 77) will 
consume major amount of the electron energy. 

In Fig. Q] the non-equilibrium spin current as a function of 
the amplitude of the external driving is shown for different 
values of the z-projection of the magnetic field. For small val- 
ues of the driving amplitude and small magnetic fields it is 
seen that if the magnetic field has the same direction as the 
z-axis, the spin current decreases, while the opposite direc- 
tion of the magnetic field amplifies the spin current. This be- 
havior can be physically understood from Eq. ( T54b . Positive 
values of Hq can be equivalently considered as larger values 
of fajjQ, that is of the distance between the transverse modes. 
This in turn leads to a decrease of the transition probabilities 
which suppresses the spin current. On the contrary, negative 
values of Ho correspond to smaller values of ftwo leading to 
an increase of the transition rates and thus the spin current is 
enhanced. Another physical explanation is that the magnetic 
field aligns the spins along its direction. Therefore, when Ho 
is positive or negative the spins are forced to point in the di- 
rection of the z-axis or in the opposite direction, respectively. 
The spin current gets more polarized in the direction of the z- 
axis for Ho > or in the opposite direction for Ho < 0. As a 
consequence its magnitude decreases for Ho > or increases 
for Hq < since it was polarized in the direction opposite to 
the one of the z-axis in the absence of the magnetic field. 

The same dependence of the spin current on the magnetic 
field with small values of its magnitude (as well as for a small 
value of the driving amplitude FL = l.Ofiwrj) is found in 
Fig. |2] in view of its dependence on the spin-orbit interaction 
strength fe so . Again for Ho > the magnitude of the spin- 
current gets smaller and for Ho < it gets larger. Addition- 
ally, one can see that the presence of the magnetic field does 
not change the locations of minima and maxima of the spin 
current as a function of fc so . This has the following physical 
explanation. The minima and maxima in Fig. |2] are related 
to the periodicity of the energy spectrum in the k-space. In 
terms of the band energy versus the quasi-momentum k de- 
pendence RSOI produces a horizontal (that is the energy of 
the bands does not change) split of the energy bands as well 
as their hybridization. Due to the periodicity this split can be 
minimal or maximal for some values of k so which leads to the 
corresponding minima and maxima in Fig. [2] The role of the 
hybridization is that the split is never zero and thus the minima 
of the spin current are not exactly equal to zero. In contrast to 
this horizontal split the magnetic field produces a vertical (that 
is along the energy axis) split and it also produces hybridiza- 
tion. This vertical split is not correlated with the periodicity 
of the energy bands in the k-space and thus the locations of 
minima and maxima remain untouched by the magnetic field. 

However, the picture explained above is only valid for small 
values of the driving amplitude F and magnitude of the mag- 
netic field I Ho I where the spin current has a linear response to 
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Magnetic field, g^ B H Q [ h co Q ] 

FIG. 3: (Color online) Non-equilibrium spin current, J^L e> g, as a 
function of the magnetic field, gfisHo, for different values of the 
amplitude of the driving force, F. The other parameters are as in 
Fig.ffl 



the magnetic field. When \Hq \ increases further, the spin cur- 
rent depends non-linearly on Ho and a complicated interplay 
between the magnetic field, driving and dissipative processes 
develops. This dependence of the spin current on the magnetic 
field is depicted in Fig.[3]for different values of the amplitude 
of the driving force. In order to stay within the validity of our 
model, where only the first two transverse modes are opened, 
the magnitude of the magnetic field must satisfy the condition: 



gHB\H \ < 0.5(ftwn + Ae 7;4 ,i), (61) 



where Ae 7;4 ,i = £ fcB [Ci ; i(fe) - 4°0;i( fc B)] /N. For the 
values of the parameters used to obtain the numerical results 
we have Ae 7; 4.i = — 0.07ftwo- Thus gfi-s\Ho\ ^ 0.465fto;o. 
As it can be seen from Fig.[3]the magnitude of the spin current 
decays for large positive values of Hq- This happens because 
the distance between the Bloch sub-bands becomes large and 
thus the transition processes are less probable. For a certain 
negative value of Ho the magnitude of the spin current has 
a maximum after which it starts to decrease and vanishes at 
some point H^ < 0. After this point and for Ho < Hq 
the spin current reverses its sign and its magnitude increases 
again. This behavior clearly demonstrates that the magnetic 
field can, without changing its direction, act in phase (i.e. de- 
stroy the spin transport) with the dissipative processes as well 
as out-of-phase (i.e. intensify the spin kinetics) with them. 
Mathematically it comes from the fact that in Eq. d32b for the 
transition rates the magnetic field Ho and the imaginary part 
of the twice integrated bath correlation function Im[Q(r)] en- 
ter the arguments of the same trigonometric functions. This 
is clarified by Eq. d37l i appropriately rewritten below for the 
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FIG. 4: (Color online) Non-equilibrium spin current, J^L e s, as a 
function of the magnetic field, gfisHo, for different values of the 
viscosity coefficient, r\. The driving amplitude is F — 1.0hojo/L. 
The other parameters are as in Fig.Q] 



case £' = 1, £ = 4: 



X COS 



X J 



ft 2 

/ Ag 7;4 ,i , , 2gfi B H Q \ L 2 

\^^ + " + ^ ( — J T -T gi(r) 

2FL . /fir 



x (62) 



sin 

nn V 2 



where Q R (r) ee Re[Q(r)], Qi(r) ee Im[Q(r)]. The phys- 
ical explanation of why in our system the magnetic field in- 
teracts only with the friction part of the dissipation and not 
with its noise part is rooted in the roles which the magnetic 
field and dissipation play for quantum coherence. On the one 
side quantum coherence in a dissipative system dies out due to 
the noise part of the Feynman- Vernon influence weight func- 
tional. On the other side, within the Feynman path integral 
formalism, we see that in our system a transverse in-plane uni- 
form stationary magnetic field cannot produce the additional 
phase due to the integral of the vector potential along the 
Feynman paths (see Appendix let. Thus in our system quan- 
tum coherence is totally insensitive to the magnetic field and 
as a result cannot interact with the noise part of the Feynman- 
Vernon influence weight functional. 

The mutual impact of the magnetic field and quantum dis- 
sipative processes on the spin current in the system is shown 
in Fig. 2] where the spin current is plotted versus the magnetic 
field, Ho, and the viscosity coefficient, 77, plays a role of a pa- 
rameter. Again for large positive values of Hq the spin current 
vanishes. As expected, the spin current gets smaller if the dis- 
sipation in the system gets stronger. When the dissipation gets 
weaker (77 = 0.05 and 77 = 0.03 curves) the oscillations of the 
spin current become observable. These oscillations are related 
to the interaction between the magnetic field and driving and 
can be described in terms of the photon emission/absorption 
processes 36 since changing Hq is equivalent to changing the 
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Magnetic Field 

FIG. 5: (Color online) Contour plot of the non-equilibrium spin 
current, J^_. e ,s [Luio], as a function of the magnetic field, gfisHo 
[hwo], and viscosity coefficient, 77. The other parameters are as in 
Fig.E 



distance between the corresponding Bloch sub-bands. 

The minima in Figs.[3]and|4]at negative values of Hq appear 
as a result of a cooperative action of the orbit-orbit coupling, 
confinement, magnetic field, driving and dissipation. Its loca- 
tion changes when the strength of the driving and dissipation 
varies. 

For completeness in Fig. [5] we also show the spin current 
as a contour plot using the variables Ho and 77. The main 
effect of the interaction between the electrons and external 
environment is the electron dressing. The dressed electrons 
are heavier and as a result less mobile. Since the spin degree of 
freedom is carried by these dressed electrons, the spin current 
decreases when the viscosity coefficient grows. 



VIII. CONCLUSION 

In conclusion, we have studied averaged stationary quan- 
tum transport in a driven dissipative periodic quasi-one- 
dimensional (ID) system with Rashba spin-orbit interaction 
(RSOI) and placed in a transverse in-plane uniform stationary 
magnetic field. For the case of moderate-to-strong dissipation 
it has been shown that the averaged stationary charge trans- 
port is well suppressed as soon as it is restricted within the 
Bloch sub-bands grown out of the same Bloch band of the cor- 
responding truly ID problem without RSOI. However in the 
same situation the averaged stationary spin transport is acti- 
vated. The analytical expression for the spin current has been 
derived and its behavior as a function of the driving param- 
eters, dissipation, spin-orbit interaction strength, orbit-orbit 
coupling strength and a transverse in-plane uniform station- 
ary magnetic field has been analyzed. Our results on the spin 
transport in the system have been presented and thoroughly 
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discussed. It has been found that the spin current as a function 
of the magnetic field shows a highly non-trivial dependence 
for different values of the dissipation and driving parameters. 
In particular, increasing the magnitude of the magnetic field 
does not always lead to a monotonous response in the magni- 
tude of the spin current. The magnitude of the spin current can 
have maxima after which its dependence on the magnitude of 
the magnetic field changes to the opposite one. Moreover, the 
spin current as a function of the amplitude of an external lon- 
gitudinal ac electric field has reversals of its direction when 
the system is placed in a finite transverse in-plane uniform 
stationary magnetic field. Also as a function of this magnetic 
field the spin current changes its direction at finite values of 
the amplitude of the ac electric field. Such behavior is un- 
doubtedly related to a deep correlation between the dissipative 
processes and magnetic field effects in the system. 
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APPENDIX A: EIGEN- VALUE STRUCTURE OF THE 
COORDINATE OPERATOR IN A SUBSPACE GENERATED 
BY BLOCH STATES OF A FINITE NUMBER OF BANDS 

In this appendix we consider a physical property with the 
corresponding quantum mechanical operator which when op- 
erating in the Hilbert space has a continuum spectrum and 
show how this continuum spectrum can turn into a discrete 
one under a certain restriction of the Hilbert space. To be spe- 
cific we constrict the Hilbert space to a subspace using some 
of the Bloch states and consider how the coordinate operator 
changes its spectrum. 



1. Introduction 

In many problems of condensed matter theory one is not 
usually interested in the full band structure of a solid but rather 
in a few bands most important for the relevant physics of a 
system. For example in metals one or a few bands with an 
energy range containing the Fermi energy are most important 
since the main contributions to transport properties come al- 
most only from those bands. Taking into account only a few 
Bloch bands leads to a restriction of the Hilbert space to a 
subspace which is then used to describe physical properties. 



2. Truncation of the Hilbert space using Bloch states 

Let H be the Hilbert space of all possible states and let us 
choose in this space the basis of Bloch's states {\l, fee}}: 

(x\l,k B ) = e ikBX UL kB (x), 

ui, kB (x + L) = ui t k B (x), (Al) 

V k B eB.Z., 1 = 1,2,... , 

where L is the period of the Bloch amplitude ui : k B {x) and 
B.Z. stands for the first Brillouin zone. 

Any vector \ip) £ H represents a linear combination 

oc 

H>)=Y, E c l,kn\hk B ). (A2) 

1=1 fc B GB.Z. 

Another basis \a) is obtained using a transformation 

\a) = U~ 1 \l,k B ), Vk B eB.Z., 1 = 1,2,... , (A3) 

where U is an arbitrary unitary operator. 

Let us consider an operator O corresponding to an observ- 
able O. Its matrix representations in the two bases ( IAU and 
dA3l are 

Ob = (l',k B \6\l,k B ), 

V k Bl k' B e B.Z., 1,1' = 1,2, ... , (A4) 
O a = (a'\G\a), Va,a'. 

The eigen-values {Xi} of the two matrices dA4l ) are the same 
and represent all possible values of the observable O. 

Now let us consider a subspace S C TL generated by 
Bloch's states corresponding to a finite number, A^b, of bands. 
A vector \ip s ) € S has the form: 

|V>>=2 Yl c kte\k,k B ). (A5) 

i=l fc B £B.Z. 

In this subspace the operator O has the matrix representation: 

oi = {k>,k B \6\h,k B ), 

V k B 6 B.Z., %,%' = 1,2,. . . N B . 

Now the eigen-values {A^} of JA61 > do not represent all possi- 
ble values of the observable O but they only give approximate 
values of some of them. If the operator O corresponds to a 
continuous observable with the spectrum from — oo to oo, the 
eigen-values {Af } are some of the eigen-values {A;}, that is 
in this case {A^} C {A;}. 

A new basis {a } of the subspace S is related to the Bloch 
one as: 

\a s )=Ui%,k B ), 
Vfc B eB.Z., i = l,2,...,N B , 

where now Ug 1 is not an arbitrary unitary operator, but a uni- 
tary operator with the following property: 

U s ■■ \v) eS U s \v) eS, V \v) g S. (A8) 
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In this case the matrix 

Of = (a' s \6\a s ), \/a' s ,a s 



(A9) 



has the same set of eigen-values {A„ } as the matrix O b in 
dA6t . 



3. Example: coordinate 

Let us specify the observable O from the preceding section 
to be particle's coordinate q with the corresponding operator 
denoted as q. We consider the operator q in the subspace S. 
Its matrix with respect to the Bloch basis is 



<7b = (h',k' B \q\h,k B ), 

V k B ,k' B e B.Z., = 1,2, 



(A10) 



Let us choose the translational operator as the unitary oper- 
ator Us from the preceding section, that is 



U s (a) 



(All) 



It is obvious that for an arbitrary value of a the operator Us (a) 
does not satisfy the property dA8l ). However, in the case a = L 
a Bloch state \l, k B ) is translated into a Bloch state with the 
same I, k B and thus (IA8b is fulfilled. Hence, the matrix 



(A12) 



qi = (l i ,,k B \Us(L)qUs 1 (L)\l i ,k B ), 
Vk B ,k' B 6B.Z. ) i,i / = l,2,...,JV B 



has the same eigen-values as the matrix q B in ( lAlOt . But due 
to the equality 



U s (L)qU s \L) = q + L 
the two matrices q B and q B are related as follows: 

(l^k'^Usi^qU^im^kB) = 

= (k',k B \q\h, k B ) + LSi'^Sk: - feB , 



(A13) 



(A14) 



V k B ,k B G B.Z., = 1,2,. 



,N B 



From ( IA14b it follows that the eigen-values of the matrix q B 
are invariant under a shift equal to jL with j being an integer. 
That is for any Xf G {Af } there exists G {A„} such that 



JL + K 



(A15) 



Let us denote through {df } those eigen-values of the dis- 
tance between which is less than L, 



\d$-dp\<L, V df,d°, e{d%}, 



(A16) 



and which are in the zeroth elementary cell. Then each of 
the eigen-values {Af } of the matrix q B is obtained from its 
corresponding eigen- value d m G {(if } by a shift jX with a 
proper integer j. It means that each elementary cell contains 
the same number of eigen-values of the coordinate operator. 
Since the total number of the eigen-values {Af } is equal to 



N B N where N is the number of the elementary cells, we con- 
clude that there are N B eigen-values in each elementary cell. 
This gives us the final expression for the eigen-values of the 
matrix q B {N is even to be definite): 



A£j =jL + di,m = l,2,...,N B , 
N N 

J = -y + l,-y + 2,...,-l,0,l,. 

and TV — > oo afterwards. 



JV N (A17) 
2 ' 2 



APPENDIX B: POLARIZATION OF THE SPIN CURRENT 
ALONG THE x AND y DIRECTIONS. CONVENTIONAL 
SPIN CURRENT 

In Ref. [34] it was shown that in an isolated system (without 
any external magnetic field) the only non-vanishing spin po- 
larization is along the confinement direction, that is along the 
z-axis. A natural question is then what is going on in an open 
driven system in a uniform stationary magnetic field applied 
along the z-axis. The external force © and bath Hamiltonian 
(|6]l couple to the longitudinal orbital degree of freedom, that 
is to the x-coordinate of our system. Because of the spin-orbit 
coupling the external force and bath affect the spin dynam- 
ics of electrons in the quasi- ID system. The magnetic field 
also influences the spin transport. Can it then happen that the 
longitudinal spin current acquires components polarized along 
the x and y axes? Below we show that the components Jg ' v (t) 
of the spin current ( fTTT > identically vanish. 



1. Longitudinal spin current components Jg y (t) polarized 
along the x and y axes 



The expressions for the spin currents 



Js' V (t) = ^{a x . y xp{t)) (Bl) 

can easily be found using the ct-DVR basis {\(, m, j, f} 7J } 
introduced in SectionUm 



x Be( y ,j{C,m,j,a' = +l\W(t)\(,m,j,a = -l) 7 ,i), 



(B2) 



Jf (f) = -2-Ti-b J2 ( mL + d 



x Im( 7J -(C, m, j,a' = +l\W(t)\C,m,j,a= -l) 7 ,j), 

(B3) 

where we have explicitly written the trace over the bath de- 
grees of freedom in order to work further with the cj-DVR 
matrix elements of the full statistical operator W(t). 
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2. Selection rules for the cr-DVR matrix elements of the full 
statistical operator 

It turns out that the case of a harmonic confinement allows 
one to formulate selection rules for the cr-DVR matrix ele- 
ments of the full statistical operator. These selection rules are 
very useful for understanding some of the properties of the 
spin transport. 

To find the selection rules mentioned above let us decom- 
pose the Hamiltonian H in (O into 



H = H n 



where 



H = 



h 2 k 2 

2m 



h 2 h. 



h 2 K. 



m 



U(x)[l + 1 — 



{o x k z - a- z k x )-g^Ba-zH Q = 
■ {&xkz — &zk x ) i 



(B4) 



(B5) 



(B6) 



and k' x = k x — gfiBHom/h 2 k so . The full statistical operator 
has the form W(t) = U(t, t )W(t )& (t, t ), where the evo- 
lution operator U(t, to) is given as the time-ordered exponent 



U(t,t a ) = Texp 



dt'H full (t') 



to 



°° / [\ n r t ^ 
= J2 / dtiH{uu(t n ) ■ ■ ■ H mi (h). 



n=0 



operator): 

H{ u ii(t3)H{ u u(t2)Hf u ii(ti) = Hr_z + H^_ Z (H + 

+ Hext(h) + ffbath)+#R-z(#0 + Hext(t 2 )+ 

+ -ffbath) #R-Z + Hr^z(Hq + ffext(*2) + -ffbath) X 

X (H + fl«t(*l) + ffbath) + (ffo + H cxt (t 3 ) + 

+ -ffbath) H-R-Z + {Ho + H ext (t3) + Hhat}i)Hji-z'X (B8) 

x (Ho + H ext (ti) + H hath ) + (H Q + H cxt (t 3 ) + 

+ -H r bath)(-ffo + H cx t(t 2 ) + 7?bath)-ffR-Z + (Hq + 
+ H cxt (t 3 ) +-ffbath) (Ho + H ex t(t2) + -ffbath) (#0 + 
+ H cxt (ti) + -ffbath)- 

Since for a harmonic confinement all the factors (Ho + 
H ex t(tk) + -Hbath) and k x couple states with indices j and 
j' only of identical parity and the factors k 2m+1 couple states 
with indices j and j' only of opposite parity, we conclude that 
the matrix elements 7 j (£, m, j, a' = +l\W(t)\£,m,j ) a = 
— l) 7 .j, being diagonal in j, are equal to zero: 

yJ ((,m,j,a' =+l\W(t)\(,m,j,a= -l) 7j = 0. (B9) 



The selection rules jB91 > represent a specific property of sys- 
tems with a harmonic confinement. From jB9i one gets 



■#"(*) = 0. 



(BIO) 



In spite of the fact that this result is only valid for the case of a 
harmonic confinement it is still general in two respects: 1) it is 
valid not only for the stationary state but for all times t to; 
2) the external force F(t) is arbitrary. 



3. Role of the spin current definition 



Only the terms of Hr_z with odd powers contain the spin 
operators. These terms are linear in a x and a z or bilinear 
in these spin operators which is equivalent to being linear in 
a y . Contributions to the matrix elements TjJ - (£, to, j, a 1 — 

+l\W(t)\(,m,j,a = — l) 7 ,j come from the first order terms 
in a x . It is easy to see that these terms represent products 
of the factors (Ho + H ext (tk) + Hbath) ordered chronolog- 
ically (we mean the chronological ordering on the Keldysh 
contour 44 and thus operators from U'(t, to) are also included 
under this terminology), an odd number of factors k z dis- 
tributed in between (Hq + H ext (tk) + ifbath) in all possi- 
ble ways and a number (even or odd) of factors k' x also dis- 
tributed in between (Ho + H ex t(tk) + -ffbath) in all possible 
ways. Such a structure is related to the fact that the Rashba- 
Zeeman Hamiltonian, £/r-z, is bilinear in the operators a x 
and k z - To clarify our above statement we write down the 
third order term coming for example from U(t, to) (a similar 
result is obtained for products which are composed from dif- 
ferent, U(t, t ), 17* (i, t ) or W(t ), parts of the full statistical 



In light of the mathematical formalism of the this appendix 
it is now convenient to discuss the difference between the con- 
ventional spin current definition and the definition of the spin 
current used in our work, that is the definition introduced by 
Shi et al£^. We will consider the z-polarized components of 
the spin currents obtained from the two definitions. The con- 
ventional spin current operator and the conventional spin cur- 
rent will be denoted as Jg onv (t) and Jg onv (i), respectively. 
The spin current operator and the spin current which are used 
in our work will be denoted as Js(t) and Js(t), respectively. 

The two definitions and the difference between them are 



J s (t) = j f (a z x), J^(t)=a z ^, 
J s (t)-J^(t) = ^-x. 



(Bll) 



One easily finds that 



dt 



2hk. 



so ~ i 

GyK Z . 



(B12) 
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Thus the relation between the spin currents is 

4 onv (t) = Js(t)+ 



2hk 



i — Tre > (mL + cL- 



C,m,j 



x (7,j(C> m >i' cr ' = + l \k z W(t)\C,m,j,a = 

- i,j(C>™>,j,a' = -l\k z W(t)\(,m,j,<T = +l) T ,j). 



(B13) 



The second term in Eq. ( IB 13b can be finite for our system. To 
show this we consider the product i?f u ii(i3)-ff full (^-ff full (ii) 
in Eq. ( IB8b . This product contains for example the term H s 
where H is given by Eq. ( IB4b . We can write this term as 



2u \ 2 



h 2 k 



h 2 k s . 



H 3 = i? n 3 + ( — ^ I Hnk 2 - ^ [ff 2 (**fc 



m 



a z k' x ) + H {a x k z - a z k x )H ] - 



h 2 h. 



\o x k z 



a z k' x )Hl 



h 2 k. ^ 3 



(a x k z - v z k' x )k 2 



h 2 k ^ 2 



k' x H k z )+k 2 HP 



[k z H k z + k' x H k' x + ia y (k z H k' x 



(B14) 



From Eq. ( IB 14b we see that the operator k z W(t) has terms 
like 



k z a y (k z H k' x - k' x H k z ), (B15) 



which are even with respect to k z and odd with respect to a v . 
Therefore, in general we have 

■yj(C,m,j,cr' = +l\k z W(t)\(,m,j,(j = -1) 7)J — 

- -t,j(C,m,j,(r' = -l\k x W(t)\C,m,j,a = +l) 7 j ^ 0, 

(B16) 

which means that the two spin current definitions are different 
in our problem. The physical reason for this can be understood 
from Eq. ( |B15l l. The term given by Eq. dB15b is finite since 
k x and k z do not commute with Hq. It happens because of the 
presence of both the periodic potential and the confinement as 
it is obvious from Eq. ( IB5b . Thus we conclude that unlike free 
Rashba electrons the two definitions of the spin current are 
not equivalent for Rashba-Bloch electrons with a transverse 
confinement. 

As one can see from Eqs. (fT3l . (f32l > and (l43l in the in- 
sulating limit Js(t) — > 0. This is just a consequence of the 
fact that the spin current definition which we use represents a 
full derivative. It is quite reasonable from the physical point of 
view that the spin ratchet effect being a transport phenomenon 
is absent in insulators. However, the conventional definition 
of the spin current is not a full derivative. The spin current 
Jg onv (i) differs from the spin current Js(t) by the second 



term in Eq. (IB 1 31 >. There is not any general physical reason 
for this term, averaged over one driving period, to vanish in 
the insulating limit at long times because it is not proportional 
to the time derivative of the averaged populations of the states 
but it is proportional to the averaged non-diagonal (in the spin 
and transverse mode subspaces) elements of the reduced sta- 
tistical operator. These averaged non-diagonal elements can in 
general be finite in insulators. The spin ratchet effect obtained 
from the conventional definition of the spin current could then 
take place in insulators which to our opinion would be un- 
physical. 



APPENDIX C: EIGEN-ENERGIES AND EIGEN-SPINORS IN 

THE PRESENCE OF ORBIT-ORBIT COUPLING AND A 
UNIFORM STATIONARY MAGNETIC FIELD ALONG THE 
z-AXIS 

In Ref. 34 periodic structures formed in a 2DEG with RSOI 
have been considered. However the influence of an external 
homogeneous stationary magnetic field on the energy spec- 
trum has not been studied. Here we generalize the results of 
Ref. l34j to the case of a uniform stationary magnetic field ap- 
plied along the z-axis. Afterwards we discuss the orbit-orbit 
coupling introduced in the main text in Eq. (0. 



1. System's Hamiltonian and the general eigen-value equation 

For an arbitrary potential V(z) (not necessarily confine- 
ment) along the z-axis and a uniform stationary magnetic field 
applied along the z-axis (2DEG is in the x — z plane) the 
Hamiltonian reads 



h 2 k 2 

2m 



H 



V(z) 



h 2 K 



(Cl) 



■U(x) - g^j, B o- z H . 



In Eq. ( ICU Hq is the z-component of the magnetic field H = 
(0, 0, Hq) and the Landau gauge, A = (—Hoy, 0,0), has been 
chosen. Additionally we have used the fact that in a 2DEG 
y = 0. This choice effectively gives only the Zeeman term. 
The eigen-states of Hamiltonian ( IC 1 b are Bloch spinors with 
the spinorial amplitude given as (see Ref.[34h 



ui,k B ,r,(x;j,a) = UL kB+akso (x)9i tkB ^(j,a), 



(C2) 



where ui y k B (x) is the Bloch amplitude of the correspond- 
ing truly ID problem without the magnetic field and without 
RSOI, and 0i,k B ,ri(j, a ) is the eigen-spinor. This eigen-spinor 
is obtained from the solution of the eigen-value equation for 
Hamiltonian ( ICU : 



h 2 kl 



2m 



h 2 k 1 (C3) 



£l,r,(kB)0l,k B ,r,(j,<j). 
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2. Harmonic confinement 

The case of a harmonic confinement is characterized by the 
following matrix elements of the operator k z ; 



(j\k z \j') = ±iS : 



(J 



2 "F 2, 



majQ 



2Ii 



and eigen-energies ef : 



e) = fuv [ j 



(C4) 



(C5) 



Therefore the only change in comparison with Ref . 34 is in the 
diagonal matrix elements of the Hamiltonian. Reproducing 
the same calculations as in Ref. 34, that is taking into account 
only the first two transverse modes (j — 0, 1, a = ±1, rj = 
1,2,3, 4), one finds that the only change in the final results for 
the eigen-energies and eigen-spinors consists in replacing the 
function e ; _ (k B ) with 

_,, „ , £/ (0) Ob + fc so ) - e| 0) (fcB - fcso) 
£ ; (k B ;H ) - - 

where we have explicitly shown the dependence on the z- 
component H of the magnetic field. The expressions for the 
eigen-energies and eigen-spinors written through the function 
ef(k B ) in Ref. are unchanged. Also the structure (that is 
the zero and non-zero components) of the four dimensional 
eigen-spinors is the same. 

The time reversal symmetry is now broken and as a result 
the symmetry relations between the eigen-energies and eigen- 
spinors hold only if one simultaneously changes the direction 
of the magnetic field. For the eigen-energies we have: 



£i,r)=l(fcB; H ) — £i,r)=2( — k B ; ~Hq), 



£i^=3(k B ; H ) — £i^=i(— k B ; —Ho)- 
For the eigen-spinors the symmetry relations are written as: 



(C7) 



Kk B ,v=iU = {0, 1}, o = {+1, -1}; ff ) = 
= 0i,k B ,v=2tt = {0, 1}, a = {-1, +1}; -H ), 
%k B , v = 3 (j = {0, = {-1, +1};H ) = 



(C8) 



where it is also taken into account that the z-projection of the 
spin operator (and as a result its eigen-values) changes its sign 
under the time reversal. The only non-vanishing polarization 
is again the one along the confinement (and also magnetic 
field) direction. The symmetry relations for its components 
are: 



Pi%i A (kB-,H ) 



-p, 



i,77=2,3 



(-ha; -Ho) 



(C9) 



Finally, we would like to note that since for the model with 
the first two transverse modes an operator even with respect 
to z is effectively diagonal, the results obtained above remain 



valid with the following change. The corresponding truly ID 
problem without RSOI and transverse confinement has to be 
solved now not for the periodic potential U(x) but for the 
periodic potential U^j(x) = U(x)[l+jh(j + l/2)/mwo-L 2 ]. 
Thus the solution of that truly ID problem acquires a depen- 
dence on the transverse mode quantum number j through 
the periodic potential dependence on that quantum number: 



(fe) 



This does not 



change the structure (location of zero and non-zero entries) 
of the resulting 4x4 matrix which is thus diagonalized in 
the same manner as in Ref. |34j. We label the eigen-energies 
and eigen-spinors obtained from this diagonalization as 
£ T ,i,ri{k B ]Ho) and $ j; i : f. ntri (j,a;Ho) to stress their depen- 
dence on the orbit-orbit coupling strength 7. The symmetry 
relations (IC7b - dC9l > are, of course, unchanged. 



APPENDIX D: BLOCH STATES IN THE DVR 
REPRESENTATION 

The scalar products ~j (£, m\l, k B ) 1 .j are nothing else than 
the Bloch states of the corresponding truly ID problem with- 
out the magnetic field and without RSOI in the representa- 
tion of the coordinate operator x operating on the subspace 
S C H, see Appendix lAl Thus using the eigen-values (given 
in Appendix[A} of this coordinate operator we have: 

7 j(C,m|l,A? B ) 7 j = e i *»( mi+ ^«)«^ fcB (d y . CJ ), (Dl) 

where we denoted the Bloch amplitude with the abbreviation 
DVR in order to stress that it originates from the discrete vari- 
able representation and differs from the one which originates 
from the continuum variable representation. 

The difference of the squares of the absolute values of the 
hopping matrix elements, |A^* e J (m)| 2 and lA™^' (m)| 2 , 
in (l52l can now be expressed in terms of the DVR Bloch am- 
plitudes as 



|A;^(m)| 2 -|A;f 4 CT /(m)| 



0^0 \ ^ 
1 ^ 

k B ,k' B 



sin[(fc B - k B )L]lm[F TtkE 



(D2) 



where we have introduced a function F~. 



defined as 



^7:fcB,fe^- M 7,0;l,feB + feso( rf 7;l:0)" 7 4:l,fc^_fc 3o (rf7;l-l) X 

x ra, fcB - fcao (^;i,i)<^, fc ' +fcso (^;i,o)]*. 



(D3) 



The function F^.^^y has two useful properties which di- 
rectly follow from its definition dD3l ). The first property comes 
from the fact that F 7 -k B .k' B is rea l if the Bloch amplitudes are 
real: 



Im[<y£ fcB (d 7;1J )]=G 



0. (D4) 



The second property is that F 1=0 .k B .k B is an even function 
in both of its arguments. Indeed, when 7 = 0, we have 
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u i,j;i,kB ( x ) = u Lk B (x), d T ^j = d^, that is property one gets 

x [u^Jd^u^Jd,)]*. iD " Im[F^ 0;kB ,_ k , B ] = Im[F 7=0;fcB , fe , ], 

One then finds from Eq. dD5l l that F 7=0 ;-fc B .fc B = -F7=0:fc B ,fc B which means that Im[F 7=0 ;fc B ,fcJ is even in fc B and fc^. The 

and F 7= o ; ft B) _fc' = Pf=o-,k B ,k' ■ As a consequence, from this same is also valid for Rc[F 7=0 :/c B ,fc' ]■ 
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